Efficient implementation of adaptive P1-FEM in Matlab
نویسندگان
چکیده
We provide a Matlab package p1afem for an adaptive P1-finite element method (AFEM). This includes functions for the assembly of the data, different error estimators, and an indicator-based adaptive mesh-refining algorithm. Throughout, the focus is on an efficient realization by use of Matlab built-in functions and vectorization. Numerical experiments underline the efficiency of the code which is observed to be of almost linear complexity with respect to the runtime. Although the scope of this paper is on AFEM, the general ideas can be understood as a guideline for writing efficient Matlab code. 2010 Mathematical subject classification: 68N15; 65N30; 65M60.
منابع مشابه
High performance implementation and conceptual develop- ment of the wave based method for steady-state dynamic analysis of acoustic problems
The wave based method (WBM) is considered for the steady-state dynamic analysis of acoustic problems as an alternative for the finite element method (FEM). The WBM exhibits better convergence properties with respect to model sizes, but the computation (CPU) times obtained by the current MATLAB implementation are large compared to commercially available FEM codes. Recently, a high performance FO...
متن کاملModeling Ionic Polymer-Metal Composites with Space-Time Adaptive Multimesh hp-FEM
We are concerned with a model of ionic polymer-metal composite (IPMC) materials that consists of a coupled system of the Poisson and Nernst-Planck equations, discretized by means of the finite element method (FEM). We show that due to the transient character of the problem it is efficient to use adaptive algorithms that are capable of changing the mesh dynamically in time. We also show that due...
متن کاملVectorized Matlab codes for linear two-dimensional elasticity
A vectorized Matlab implementation for the P1 finite element is provided for the twodimensional linear elasticity with mixed boundary conditions. Vectorization means that there is no loop over triangles. Numerical experiments show that our implementation is more efficient than the standard implementation with a loop over all triangles.
متن کاملImplementation of $hp$-adaptive discontinuous finite element methods in Dune-Fem
In this paper we describe generic algorithms and data structures for the implementation of hp-adaptive discontinuous finite element methods in the Dune-Fem library. Special attention is given to the often tedious and error-prone task of transferring user data during adaptation. Simultaneously, we generalize the approach to the restriction and prolongation of data currently implemented in Dune-F...
متن کاملDUNE-FEM. A general purpose discretization toolbox for parallel and adaptive scientific computing
DUNE-FEM is a free discretization toolbox for parallel and adaptive scientific computing based on DUNE. The implementation of discretization schemes such as finite elements, finite volumes or discontinuous Galerkin schemes is based on abstractions that are very close to the mathematical description of the underlying methods. In this contribution we will give a compact overview on the design and...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Comput. Meth. in Appl. Math.
دوره 11 شماره
صفحات -
تاریخ انتشار 2011